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Abstract 



The density of states for an extended MIT bag model is studied numerically 
by using a parameterized smooth representation which provides the best fit to 
the numerical data. It is found that the mass dependence of the surface term 
in the density of states agrees with that derived from multi-reflection theory 
calculation. The mass dependence of the curvature term in the density of 
states is extracted for finite values of the quark mass. The scaling properties 
of the data, which is assumed in the parameterization, are studied. The 
difference between the parameterized smooth representation determined by a 
best fit and the results derived from direct smoothing of the numerical data 
are discussed. The fiuctuations in the density of states are also discussed. We 
provide a smooth representation of the density of states of a spherical cavity 
for non-relativistic particles. 

PACS number(s): 12.39.Ba, 12.39.Ki, 21.10.Ma and 12.39.x. 



I. INTRODUCTION 

Since the bag model suggested by the MIT group in the 1970s, many other bag models 
including the chiral bag [Q, soliton bag and quark- meson coupling model which take 
the quark and meson coupling interactions into account, have been put forward with widely 
successful applications in nuclear physics. Recently, the bag model has been extended to 
include strange quarks and has been employed to study the properties of a strangelet 
0]- IP which is believed, if it exists, to be an unambiguous signature for the formation of 
quark-gluon plasma 0. 

To calculate the dynamical and/or thermodynamical quantities of spherical nuclei or a 
spherical droplet of hadronic or quark matter, in particular the strange quark matter, by 
means of the bag model, we must first look for the density of states (DOS) of the spherical 
system. The DOS of a spherical cavity in which the free particles are contained can be 
expressed as |]TU| 



where N{k) is the total number of particle states. In general, the quantum states are discrete 
and N{k) is a summation of discontinuous states. But in practical calculation, a continuous 
representation of the DOS is needed for changing the summation to an integral. It will bring 
some errors and we will discuss the problem in detail below. 

Usually, N{k) can be written in terms of the dimensionless variable kR as 

N{k) = A{kRf + B{kRf + C{kR) + Ar^(°l(fc), (2) 

where R is the radius of the bag. The first, second and third terms of the right hand side 
of Eq.(0) refer to the contributions of the volume, surface and curvature, respectively, and 
^reK^) refers to other contributions which are not proportional to (fc-R)^, {kRY or {kR). 
We will discuss the term N^^^{k) in section 3 in detail. In general, the coefficients A, B 
and C are expected to be functions of m/k, and their expressions are model dependent. For 
example, using the box renormalization of the wave function, Jaqamin et al obtained [|lT 
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A=—,B = --,C = ^. (3) 
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In their approximation. A, B and C are all constants and do not depend on the momentum. 
Of course, this treatment seems to be too rough because the plane wave is not the real 
eigenmode of the spherical bag geometry. The exact eigenf unctions of the spherical bag 
model, as shown earlier, are the spherical Bessel functions. 

Another approach which is better than the plane wave approximation is the multi- 
reflection theory (MRT) |T^, [jl3|. According to MRT, the volume term A can be fixed 
using the known total number of states for an infinite system with very large R, 

limA = A^ = |f, (4) 

where g is the total number of degeneracies. For example, g is the total number of spin and 
color degrees of freedom for a quark with flavor treated separately. The surface term B was 
approximately derived from MRT to be H 



mi \k 



where m is the mass of a quark. The curvature term C has not even been evaluated the 
present except for the two limiting cases: m — > and m ^ oo according to MRT. Madsen 
proposed [JI3[ an arbitrary form for the curvature term C that simply connects the 



above mentioned limits, namely 

~,m, 9 fl ( k m\ , k vr/c 1 

to study the strangelet. Of course, this is only an ansatz with two correct limits. 

Since the spherical bag geometry is popular for flnite nuclei or quark systems and the 
DOS plays a key role in calculations, it is of interest to study this problem by exact numerical 
calculations under such a geometry. Gilson and Jaffe |]16[ evaluated the microscopic total 



number of quark states N{k) numerically under the MIT bag boundary condition for a flxed 
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quark mass m = 150MeV. This is appropriate to the nonrelativistic hmit m > k . To our 
knowledge, the mass dependence of the DOS found by numerical calculation is still lacking. 
The purpose of this paper is to check the MRT, especially the surface term given by Eq.(^ 
and the curvature term by Eq.(|^), and to find the mass dependence of the DOS by numerical 
counting of the eigenmode under the MIT bag boundary condition. This is of value because 
our knowledge of N{k) is improved through the following analysis: 

(1) A numerical study of the mass dependence of the surface term and the curvature 
term using a smooth representation which would best fit the total number of quark states 
obtained from the numerical data. (2) Comparison of the surface term to the MRT result to 
provide a first principle's check of the theory. (3) Extraction of the unknown curvature term 
from the numerical data to provide a suitable smooth representation, instead of Eq.(^, our 
representation is not only correct at the two limits m — and m — oo, but can also be 
used in the whole mass region. (4) Study of the extent to which the smooth representation 
of N{k) can be applied. 

The organization of this paper is as follows: Section 2 contains a description of the 
procedure for numerical calculation of the total number of quark states. In section 3, the 
resulting fiuctuating numerical data are fitted by using a smooth representation for the 
same quantity to obtain the coefficients of the volume, surface and curvature terms. The 
fiuctuations and their behaviour under smoothing of the remaining part of the numerical 
data are also studied in this section. Section 4 is our conclusion, discussion and the results 
for the DOS using a non-relativistic boundary condition. 



II. NUMERICAL EVALUATION 



The quarks inside a MIT bag are free Dirac particles bouncing back and forth from the 
infinite barrier on the bag surface. Their discrete energies in a spherical bag are determined 



by the MIT boundary condition |17 



(7) 



where ip is the wave function for quarks, dV the MIT bag surface and n'^ the outward unit 
vectorQ which is normal to the MIT bag surface in its rest frame. The eigenstates of the 
quarks in such a spherical cavity can be written as 



1pnKM{r) 



(8) 



with X = qr and q the absolute value of the 3-momentum of the quark. Here, r is the 
unit vector pointing in the same direction as the coordinate vector r for the quark. The 
sub-indices n = 0,1,2, k = ±1, ±2, . . . , and M are quantum numbers besides the fiavour 
and color that completely specify a single quantum state for a quark. The radial wave 



^The convention for the metric tensor adopted here is {g^"^} = diagjl, —1, —1, — 1}. 
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functions Uuk and Vnn are proportional to spherical Bessel functions since the quarks are free 
particles inside the bag. The spin-orbital angular momentum coupled wave functions, have 
a total angular momentum j = |k| — 1/2 are y^m and y-nm- The orbital angular momentum 
/ for y^rn is j' — 1/2 if K < and j + 1/2 if k > 0. Imposing Eq.(^ on the wave function 
Eq.(|D, we get 

UnMR) = VnniqR), (9) 

which can be numerically solved for each n and n to get a discrete set of g„K. The energy e 
for a quark in such a state is 



qL + m^- (10) 

This has a one to one correspondence to g^^. It is expected that the value for each e^^ has 
to be obtained numerically for the MIT bag boundary condition. 

The approach of summation over the microscopic energy levels to obtain the thermody- 
namical potential is rather computationally expensive. For example, the thermodynamical 
potential of a strangelet can be written as 

oo oo 

n = -^ E E E (2j« + l)ln(l + e-'^(^^— (11) 

f={u,d,s} n=0 K=~co 

where (3 is the inverse temperature, the summation over k excludes k = 0, is the total 
angular momentum of the state corresponding to n, fif is the chemical potential for flavor 
/ , ric = 3 is the degeneracy due to the color degrees of freedom of quarks, and fl' is the 
contribution to the thermodynamical potential from other degrees of freedom of the system. 

In most of the cases, the quantity being summed is a smooth enough function of the 
quark energy on the scale that the DOS fluctuates. To study the gross features, we can 
introduce the concept of DOS and replace the microscopic summation by an integration 
over the quark's momentum, namely 

J2F{en.)^ / dkp{k)F[e{k)]. (12) 



Here e„K is the discrete energy level, e{k) is its continuous interpolation using Eq.(|l^), and 
p{k) is the DOS which is determined by the total number of quark states N{k). In the 
numerical calculation, the N{k) that provides the best fit to N{k) is used as a smooth 
representation of N{k). With a smooth N{k), the DOS p{k) in the above equation is then 

M.)^^. as) 



There is no guarantee that the left hand side of Eq.(|T2|) is the same as its right hand 
side in the case studied here, since the microscopic distribution of the energy levels does not 
follow a trend regular enough to avoid the fiuctuations introduced by N{k) replacing N{k). 
In fact, as is found in the following, the fiuctuations in the DOS are large and irregular. 



However, the situation can be improved by expressing Eq. ([T^) as 
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^F(e„0=/ dkp{k)F[t{k)]+<6F>, (14) 
Jo 



where < 5F > is by definition the difference between the left and right hand side of Eq.(0). 
It is easy to see that < 6F >= if 

<N-N>{k)=Y^ [N{k,) - N{h)] ^ 0, (15) 

hi 

where the summation is over the energy levels of quarks with a magnitude of 3-momentum 
in the neighbourhood of k in which F[e(A;)] does not change appreciably. The best we can 
do to find the general trend of the problem is to extract from N{k) as precisely as possible 
the smooth quantity N{k) which can be attributed to the DOS . The quantity < 6F > is 
then usually small and can be ignored. 

The total quark states with the momentum q < k for each flavor is counted according to 

oo OD 

N{k)=n,J2 E ^(A:-gn.)(2j« + l), (16) 

71=0 K= — OO 

where g„K is determined by Eq.(|^) and the step function 6'(x) = 1 if x > 0, and 9{x) = 
if X < 0. This number is, in principle, a function of three variables, namely, the cutoff 
momentum k, the mass m and the bag radius R. Since N{k) is a dimensionless number, 
it is expected to be a function of dimensionless variables too. Two such kinds of variables, 
namely, m/k and kR formed by k, m and R have been chosen in Eq.(^ while the other 
variables, for example mR, will not emerge in expressions of coefficients A, B and C because 
they are independent of R. 

The numerically counted N{k) is therefore fitted by using the following smooth trial 
functional form ansatz 

N{k) = A (^) {kRf + B (^) {kRf + C (^) {kR) , (17) 

where all coefficients A, B and C are assumed to be independent of kR. 

The detailed fitting is done in the following way. First, m and k are kept fixed (so m/k 
is a constant), and various corresponding N{k) are counted by varying R. Such a series of 
N{k) data as a function of R is fitted using the smoothed representation Eq.(p!7D to obtain 
the coefficients A, B and C for a fixed m/k. We employ "the method of minimum squares" 
to obtain the best fit. Then the value of m is changed with a fixed k and the above steps 
are repeated to obtain a different A, B and C for different m/k. 

The consistency of the procedure is checked by calculating N{k) with a series of different 
cutoff momentum k (25,50,100,200,400 and 800 MeV) with m/k and kR fixed. It is found 
that the choice of cutoff momentum will not affect the coefficients A{m/k), B{m/k) and 
C{m/k) for each m and k bin. 



III. SMOOTHED PART OF THE DOS UNDER THE BAG BOUNDARY 

CONDITION 

The results of our best fit for smooth curves by using "the method of minimum squares" 
can be summarized as follows. 
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A. The volume term 



The volume term A is found to be 

A = 7.074 X lO-'g, (18) 

which is different from 2g/97T by a difference of 10~^— 10~^. It means that A{m/k) is 
independent of m//c to a very good approximation. This result is in good agreement with 
the MRT as well as the plane wave approximation. 



B. The surface term 

The coefficient B{m/k) depends on m/k explicitly. The results of our best fitted 
B[m/k) is compared to the one from the MRT. The relative difference between B{m/k) 
and BMRri^/k) is shown in Fig. 1. It can be seen that such a relative difference is around 
1 — 2 % for most of the m/k value except at rather small m/k. Thus we come to the 
conclusion that the coefficient BMRTim/k) describes the surface correction of the DOS very 
well. 



C. The curvature term 

After we remove the volume term and the surface term , the remaining numerical data 
for N(k), which is denoted as Nj;ll{k), contains information about the coefficient of the 
curvature term C{m/k). The fluctuations of Nj:ll{k) are large. To extract the linear term, 
an average of the remaining data Nj^lKk) defined by 

, kR+SxA/2 

NTm) = ir E Nillix), (19) 

^ x=kR-5x A/2 

where 6x is the difference between the neighbouring kR points in the numerical data for 
N{k), is made first. The averaged N^ref' is fitted to a linear function by using the method of 
minimum squares to get the curvature term. The resulting coefficient C{m/k) is not a slow 
varying function of m/fc as is the one assumed in Eq.(^). The result of our best fit for A = 50 
is shown in Fig. 2 which agrees with the multi- reflection theory at the two ends, m/k = 
and m/k oo, but is a fluctuating function of m/k in between. If these fluctuations are 
smoothed out, it contains a global maximum which is different from the proposed monotonic 
function given in Eq.(^). A smooth representation of the extracted (fluctuating) function is 
the following 

C{m/k) = C{m/k) + (^) 7—^—2 (20) 



3.42(^-6.5) +100 



with C{m/k) given by Eq.(|^. This is plotted as the solid line in Fig. 2 together with the 
numerical data points. The dashed line is C{m/k). 
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D. Fluctuating term 



Most of the fluctuations N{k) reside in the remaining part of Eq.(^). It is also expected 
that the difference between the smoothing procedure of the MRT and the one adopted here 
resides in this residual term, which is obtained from the numerical data N[k) using the 
following equation 

AtW = N{k) - '^{kRf - B{mlk){kRf - C{m/k){kR), (21) 

where B{m/k) and C{m/k) are obtained from the procedure above described. We use the 
method for A^'W of Eq. to smooth N^°2 and obtain Nj^°l^{k). 

The results for A = 1, 10,50, 100 and m/k = 0.1, 100 are shown in figure 3 and figure 
4, respectively. Two trends for the numerical data are seen after subtracting the volume, 
surface and curvature terms obtained by best fitting. First the magnitude of the fluctuations 
decreases as 1/ \/A. Secondly, the magnitude of the fluctuations increases with kR. The 
fluctuations will exist for any value of A. 

This result means that iV,(o) can not in principle be treated as zero. According to Eq.(0), 
this means that the fluctuating term < SF > for any smooth function of k is not zero for the 
MIT bag model. Fortunately, it is not very large for a smooth function {N^'^l/N ^ 10~^ ~ 
10~^) and can be dropped in most of the applications involving the gross features of the 
problem. 



IV. CONCLUSION AND DISCUSSION 

The distribution of the quark states in a spherical MIT bag is calculated numerically. 
By making a best fit to the data using a smooth function of the dimensionless variables kR 
and m/k, the coefficient of this kind of power expansion is determined. It is found that the 
volume and surface terms can be represented by a smooth function of the above mentioned 
variables in agreement with the MRT. The curvature term which contains fluctuations is 
extracted from the numerical data and represented by a smooth function, as given by Eq. (pO|) . 
The fluctuations remain after smoothing, but can be ignored in most problems. 

The same treatment can also be applied to the non-relativistic case: For simplicity we 
only give the final result here. The DOS for non-relativistic particles confined in a spherical 
cavity of radius R with boundary condition 

u^{kR) = (22) 

can be found and the best smooth fit to the numerical result is 

N{k) = —{kRf - 0.1289g{kRf + l.OUAgikR) - 60.319^. (23) 

Qtt 

A non-spherical bag is needed to represent a droplet at high angular momentum, which 
is expected to be produced in relatively large number in a high energy relativistic heavy ion 
collision. This is not studied in this work. A representation and numerical study of these 
deformed bags should be carried out in the future. 
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V. FIGURE CAPTIONS 

FIG.l. The difference of the coefficient 5 as a smooth function of m/k between the 
one best fitted to the numerical counted value and the one calculated analytically. Their 
difference draw in line is almost a constant with a variation less or equals to 1% of its own 
value. 

FIG. 2. The coefficient C extracted from n'^^^^ . The solid line is a smooth representation 
of it. The dashed line represents C{m/k). 

FIG. 3. The behaviour of the fluctuating residual term Nj^'^g under the smoothing for 
m/k = 0.1. 

FIG. 4. The behaviour of the fluctuating residual term N^^^ under the smoothing for 
m/k = 100. 
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